function printout_summation_of_amplitude(m_max, CONSTS, plot_data, odn_pl)

    num_phi_points = 10051; % 50
    phi_vec = linspace(-pi, pi, num_phi_points)';
    
    % current distribution obtained numerically
    I_sum = func_I_sum(phi_vec, m_max, CONSTS); %func_I_sum(phi_vec, m_max, CONSTS);
    I_max = max(I_sum);
    ang = to_degrees(angle(I_sum));
       
    if(odn_pl)
        k0 = CONSTS.k0;
        eps = CONSTS.eps;
        eta = CONSTS.eta;
        eps_a = CONSTS.eps_a;
        a = CONSTS.a;
        d = CONSTS.d;
        c = CONSTS.c;
        Z0 = 4*pi/c;
        
        % current distribution in the case of a homogeneous magnetoplasma
        h = k0*(abs(eps*eta))^(1/4)*((1-1i)/sqrt(2)); 
        I_sum_odn = -((1i*pi*h)/(Z0*k0*log(4*a/d)))*((cos((pi-abs(phi_vec)).*h.*a))./(sin(pi*h*a))); 
        ang_odn = to_degrees(angle(I_sum_odn));
        ang_odn_correct = ang_odn;
        ang_odn_correct(ang_odn > 0) = ang_odn(ang_odn > 0) - 360;
        
        % current distribution obtained analytically
%         eps_ef = ((2i*eps_a)/((1+(eps_a)^2/(abs(eps*eta)))*(sqrt(abs(eps*eta))))) - 4;
        eps_eff = ((abs(eps*eta)+eps_a^2))/...
                    (-2i*eps_a+2*(sqrt(abs(eps*eta))));
        h1 = k0*sqrt(eps_eff)*(1-1i)/sqrt(2);
        I_sum_analytical = -((1i*pi*h1)/(Z0*k0*log(4*a/d)))*((cos((pi-abs(phi_vec)).*h1.*a))./(sin(pi*h1*a)));
        
    end
    
    if(plot_data)
        if(odn_pl)
            figure; plot(to_degrees(phi_vec), abs(I_sum)./abs(I_max), 'b-', ...
                         to_degrees(phi_vec), abs(I_sum_odn)./abs(max(I_sum_odn)), 'b--', ...
                         to_degrees(phi_vec), abs(I_sum_analytical)./abs(max(I_sum_analytical)), 'r-'); 
            title('|I_{\Sigma}/I_{Sigma, max}|'); xlabel('\phi, degrees');
            legend('plasma column numerical solution', 'homogeneous plasma', 'plasma column analytical solution');
            figure; plot(to_degrees(phi_vec), real(I_sum), 'b-', to_degrees(phi_vec), imag(I_sum), 'r');
            legend('Re(I_{\Sigma})', 'Im(I_{\Sigma})'); xlabel('\phi, degrees');
            figure; plot(to_degrees(phi_vec), ang, 'b-', ...
                         to_degrees(phi_vec), ang_odn_correct, 'b--', ...
                         to_degrees(phi_vec), to_degrees(angle(I_sum_analytical)), 'r-');
            title('angle(I_{\Sigma})'); xlabel('\phi, degrees');
            legend('plasma column numerical solution', 'homogeneous plasma', 'plasma column analytical solution');
        else
            figure; plot(to_degrees(phi_vec), abs(I_sum)./I_max); 
            title('|I_{\Sigma}/I_{Sigma, max}|'); xlabel('\phi, degrees');
            figure; plot(to_degrees(phi_vec), real(I_sum), 'b-', to_degrees(phi_vec), imag(I_sum), 'r');
            legend('Re(I_{\Sigma})', 'Im(I_{\Sigma})'); xlabel('\phi, degrees');
            figure; plot(to_degrees(phi_vec), ang);
            title('angle(I_{\Sigma})'); xlabel('\phi, degrees');
        end

    end
    
%     coef_sgs_to_si = 9e+11;
%     I_sum_0 = I_sum(5026); %func_a0(CONSTS);
%     Z = coef_sgs_to_si/I_sum_0;
%     fprintf('Z = %f + i%f\n', real(Z), imag(Z));

end